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ABSTRACT 


High voltage overhead transmission line route selection has become a vital component of the 
power transmission line project life cycle as demand for electricity has increased and global 
awareness of the environmental implications of large-scale power projects has grown. The project 
necessitates meticulous planning, engineering designs, facility installations, commissioning, and 
maintenance, as any errors made at this stage may have a flowing detrimental impact on the 
following phase, which may be too costly to fix. Using geographic information system (GIS), multi- 
criteria decision analysis approach, and environmental impact factors strategies, this study 
determines the best route from existing geographical space between Lafia 330/132/33kV Substation 
in Nasarawa State and New Apo (Pigba) 330/132/33/kV Substation in the Federal Capital 
Territory, Central Nigeria, Abuja. Land use, land cover, and other environmental data were 
obtained by performing image classification on Landsat § satellite image using ArcGIS 10.1 to 
achieve this goal. Three key terrain elements were identified to define the land cover types: built-up 
areas, vegetation, and water bodies. These layers were created to reveal the locations of significant 
settlements, geological formations, and the HV-OHTL route that avoided crossing protected areas. 
Using pair-wise comparison of analytic hierarchy process (AHP) approach, each criterion 
evaluated was assigned relative weight that shows their level of influence on the proposed power 
transmission line, and these were provided as GIS layers for visualizations. To build the cost 
surfaces assessed during route selection for HV-OHTL, weights derived from AHP method for each 
decision makers were applied to all layers to create the cost surfaces analyzed during route 
selection for HV-OHTL. 


Keywords: Geographic Information System (GIS), Multi-Criteria Decision Analysis (MCDA), 
Analytic Hierarchy Process (AHP), Least Cost Path, High Voltage Overhead Transmission Line 
(HV-OHTL), Substation 


1.0. Introduction 


The high-voltage over-head transmission line (HV-OHTL) is a dependable option for bulk power 
transmission around the world. Route selection is a difficult engineering decision making. For 
example, a balance is struck in the United Kingdom (UK) between constitutional demands, 
engineering needs, economic viability, land use land cover (LULC), and environmental stressors 
(Marshall and Baxter, 2002). This is analogous to what is required to achieve a route that is 
technically, financially, and commercially viable while also avoiding or minimizing environmental 
impact. The underlying issues in HV-OHTL route selection projects have always attracted public 
interest and sparked heated debates, with the fundamental issues frequently focusing on local effects 
such as visual configurations of associated facilities and their environmental impacts rather than the 
project's overall benefits to society. Given the scope of public and regulatory investigations, the 
environmental impact statement resulting from such plans must be impartial and transparent in the 
route selection process. The method is founded on evidence that the non-aesthetic layouts, the grade 
of visual breaks, and how environmental consequences can be minimized by well-planned route 
selection are significant issues of HV-OHTL. 
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The goal of HV-OHTL optimum routing is to efficiently transfer large amounts of electricity from the 
send end to the receiving end, where electricity is consumed at a break-even cost. For example, in the 
United Kingdom, this is frequently done at voltage levels of 132kV, 275kV, or 400kV (Scottish 
power, 2001), whereas in Nigeria, power transmission voltage levels are 330kV and 132kV. The 
influence on the right—of—way increases as the transmission voltage level rises. Transmission towers 
and lines are very hard to conceal, but there is a lot that can be done to make them less conspicuous 
and attractive. In both rural and urban areas, this includes suitable corridor site selection, landscaping, 
and screening. However, transmission line location constraints are typically greater in urban areas. 
Existing development and a lack of right-of-way may limit the best locations for HV-OHTL project. 
In rural locations, there is sufficient corridor for HV-OHTL to be built. The area enables right-of-way 
determination, project route alignment and development. 


The HV-OHTL 1s a large linear structure erected in a landscape setting, making it the most important 
in determining environmental variables that influence the best power transmission line routing. All 
those parameters that influence route decisions at varied degrees and phases of project installation and 
operations are all those parameters that exist in a transmission line right-of-way. These variables can 
have a direct or indirect impact on the path, and they might be apparent or invisible. Landscape 
environmental variables were primarily used as criteria in this study to aid in the HV-OHTL selection. 
Integrating MCDA techniques into GIS-based processes during HV-OHTL selection planning is 
helpful because it allows value judgments to be introduced into GIS-based decision-making 
procedures and feedback on their consequences for policy implementation and evaluation 
(Malczewski and Rinner, 2015). Although this combination looks to be the finest option, it is just one 
of many options available in practice (Schito et al., 2021). Electricity demand is rising internationally 
as a result of rapid population increase (International Energy Agency, 2019), and the benefits of a 
reliable and high-quality electricity supply to customers boost economic growth and social-economic 
development (Taylor, 2017). 


2.0. Methodology 


2.1. Description of study area 


The study area is located between the Federal Capital Territory (FCT) and Nasarawa State. The 
planned proposal for constructing a 330kV high voltage transmission line from Lafia's 330/132/33kV 
substation to Abuja's New APO (Pigba) 330/132/33kV substation is covered within the study area. It 
lies between Latitude 08° 20 01” N to 09° 04 00” N and Longitude 07° 16 55 E to 08° 50’ 32” E 
(Figure 1). 


2.2. Geographic Information System Software, In-Built Tools and Applications 


An advanced version of the Environmental Systems Research Institute's (ESRI) product, ArcGIS 
Desktop 10.1, was used in this investigation. It was possible to create, query, mapped, and analyse 
cell-based raster data, perform integrated raster/vector analysis, derive new information from existing 
data, query across multiple data layers, and fully integrate cell-based raster data with traditional vector 
data sources using ArcGIS Spatial Analyst. In addition, ArcGIS Spatial Analyst was integrated into 
the ArcGIS interface, allowing users to access all of the available functionalities to obtain information 
about geospatial data such as terrain analysis, spatial relationships, suitable locations, and the total 
cost of travel from one point to another (ESRI White Paper, 2001). 


Additionally, MS Excel 2014 and MS Word 2014 were employed for a variety of tasks, including 
graph production, report writing, and AHP computations. 
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Figure 1: Research study area 


2.3. Data sources and collections 


DEM, Slope, Soil, Erosion, Geology, Satellite Imagery, Land Use Land Cover, and Road Network 
layers are among the data needed for this project. 


They are made up of a variety of ages, sources, formats, qualities, and scales that have been collected 
from various organizations. The datasets and their sources and ages are presented in Table 1. 


Table 1: Datasets properties and sources 


S/N Data Names Scale /Year Data Sources Year Format 

1 Digital Elevation Model 30m Resolution 30m DEM _ downloaded from 2014 Raster 
(DEM) United states geological survey 
(USGS) (www.explore.USGS.gov) 

2 Slope 30m Resolution Generated from 30m DEM = 2014 Raster 


downloaded from United states 
geological survey (USGS) 


(www.explore.USGS.gov 


3 Soil 1km Resolution Harmonised World Soil Database 2013 Raster 
(HWSD) www.tao.org 

4 Erosion 1km Resolution K factor of Soil data (generated 2013 Raster 
from soil data) 

5 Geology 1:50,000 Nigerian Geological Survey 2005 Raster 

6 Satellite Imagery 30m Resolution Landsat 8 satellite Image 2018 Raster 
(www.explore.USGS.gov) 

7 Land Use Land cover 30m Resolution Generated from Landsat 8 satellite 2018 Raster 
Image 

8 Road network Vector data from Open street map 2018 Vector 


(www.bbbike.com ) 


2.3.1. Digital Elevation Models (DEM) 


Digital elevation models are raster-based data (e.g., a grid of squares or a Triangular Irregular 
Network, TIN) that provide a digital representation of data across a defined region at a specific scale, 
in this example 30 meters. The elevation of the route (site) varies from 65m to 1200m (Figure 2). 
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Figure 2: DEM of study area 


2.3.2. Slope (Using 30m DEM) 


Using the spatial toolbox in ArcGIS, the slope data was created from the 30m DEM (Figure 3). The 
Slope data specifies the greatest rate of change in value from one cell to the next, or a measure of 
change in surface value over distance, either in degrees or as a percentage. The flatter the terrain, the 
lower the slope value; the steeper the terrain, the greater the slope value. These two characters are 
crucial when it comes to tower design and installation. The slope classes listed in Table 2 are 
commonly recognized in engineering practice, according to (ESU, 2010). 
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Table 2: Interpretation of slope classes (ESU, 2010) 
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Figure 3: Slope classification in percentage 
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2.3.3. Soil layer 


Various soil classes were considered based on textual soil properties as it relates to the strength of soil 
present at the project site. Data for this purpose were acquired from the Harmonized World Soil Geo- 
database (HWSD version 1.2), the HWSD is composed of a raster image file and a linked attribute 
database. The raster database consists of 21,600 rows and 43,200 columns, of which 221 million grid 
cells cover the globe’s land territory (FAO, 2009). Organic carbon, pH value, water storage capacity, 
soil depth, cat-ion exchange capacity of the soil and the clay fraction, total exchangeable nutrients, 
lime and gypsum contents, sodium exchange percentage, salinity, textural class, and granulometry are 
all linked to each grid cell in the database. HWSD uses GIS software to display or query soil 
compositions in terms of user-selected soil parameters; in this case, the database was loaded into Arc 
map 10.1. The soil layer is shown in Figure 4 and it gives the major classes of soil within the region 
considering their textural properties. About 60% of the soil types were seen to be clay loam, and the 
remaining 40% consist mainly of sand and loamy soils. 
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Figure 4: Soil textural classification 


2.3.4. Geology 


Figure 5 shows the geological layer created using a Nigerian geological map. It was digitized using 
ArcGIS 10.1 software and geo-referenced for the study area. The site's geological formation reveals 
information about the material strength that should be considered during tower foundation design and 
construction. 
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Figure 5: Geological map 





2.3.5. Land Use and Land Cover 


The land cover map was created using ArcGIS 10.1 and a supervised image classification on a 
Landsat-8 satellite image using the USGS Land Use Land Cover categorization scheme (Chandra et 
al., 2014). Built-up, Vegetation, and Water Bodies were the three key land class features used to 
determine the Land Cover type. A false color composite was created using a Landsat Band 
combination of five near-infrared (SNIR), 4(Red), and 3(Green), as seen in Figure 6. Using colour 
discrimination, vegetation appears red, built-up areas seem yellow, and water bodies appear blue. A 
training data set was constructed by picking multiple known locations to represent each class of 
information, as shown in Figure 6, and then utilized to perform image classification. 


These layers display the locations of significant towns and water bodies, as well as power line route 


selection, which minimizes tower spotting in water bodies and avoids routing power lines through 
villages as much as possible. 
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2.3.6. Erosion 


For the estimation of soil erosion in the study area, parameters from the soil layer of the harmonized 
world soil geo-database (HWSD) were used, and the distribution of soil particles was taken into 
account. Although soil erosion factor K is influenced by particle size distribution, organic matter 
content, soil structure, and permeability, the silt percentage content from whole top soil content is the 
most important factor in soil erosion because silt is easily detached and tends to crust, resulting in 
high runoff rates. It has been proven that soils with a high silt concentration are the most prone to 
erosion. The erosion factor K, was calculated using (William, 1995) equation. To estimate the erosion 
factor (K) Equation] is used. 


Kysie a fesand x fei-si x forge x Shisana (1) 


Where each factors are calculated using Equation 2 — 5 





Msilt 
fesang = (0.2 + 0.3 4 See (1- T0g )p (2) 
M slit 

ee ee (3) 

fina = (— 
_ ( 0.0256. orgC Fl 
forge = orgC + exp[3.72 — 2.95. ed) 4) 

0.7.(1 — Js) 
i0 

fhisand = | 1- Ut (5) 





(1 -— 755) + exp[-5.51 + 22.9(1- 755) 


The computed values using Equations 2 — 5 are shown in Table 3. 

Where: 

m,: The percentage of sand fraction content (0.5-2 mm particle diameter) 

Mg it: The percentage of silt fraction content (0.002-0.05 mm particle diameter) 
m, : The percentage of clay fraction content (<0.002 mm particle diameter) 
orgC: The percentage of organic carbon fraction content (%) 
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Table 3: Calculation for each soil class 
Soil Textural 


Class SAND SILT CLAY Organic Content fcSand fcl_si Forge Fhisand k factor 
Sandy Loam 81 9 10 0.55 0.20 0.80 1.00 0.74 0.12 
Loam 42 36 22 1.00 0.20 0.87 0.99 1.00 0.17 
Sand 43 34 23 1.40 0.20 0.86 0.98 1.00 0.17 
Sand 73 18 9 1.40 0.20 0.89 0.98 0.92 0.16 
Sand 73 18 9 1.40 0.20 0.89 0.98 0.92 0.16 
Clay Loam 49 29 22 1.01 0.20 0.84 0.99 1.00 0.17 
Clay Loam 49 29 22 1.01 0.20 0.84 0.99 1.00 0.17 


The generated erosion map is shown in Figure 7, with average values ranging between 0.11 - 0.17. 
When erosion ratio is less than 10, it indicates non-erodible soil, where erosion ratio ranges from 12 to 
115 shows erodible soil. 
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Figure 7: Soil erosion map 


2.3.7. Road network layer 


The road network refers to the roads that exist in the research area. The proximity of the HV-OHTL 
route to the main road 1s critical because it minimizes the cost of transportation of men and materials 
during project construction, commissioning, and maintenance period. The road layer (Figure 8) was 
retrieved as a vector data type from the open street map database, and the networks were updated 
using archival satellite imagery from Google Earth Engine from 2019. The primary route in the study 
locations is depicted in Figure 8. 
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Figure 8: Major roads around study area 
3.0. Results and Discussion 


The aforementioned layers were classed as "a soil map displaying distinct soil types" and "other 
variables exhibiting different layers," respectively. In this instance, it's preferable to give the soil units 
a suitability property for evaluating different types of tower foundations. Because different soil types 
may have varied suitability’s, a classification may mix soil units of diverse types into the same 
category of soil suitability- based on how appropriate the current characteristics along the proposed 
HV-OHTL project route are. For raster calculations, the data was first transformed to usable raster 
sets. Consequently, the most suitable features for route selection were ranked as | and the worst 
features ranked as 9 in-line with Saaty scale (1980). When this method of overlaying the raster layers 
was used, the cell that has the minimum value corresponds to the most suitable cell for HV-OHTL 
project routing. Table 4 shows the assigned values to each layer. 


The reclassification values for each raster layer were redistributed using the tool in Arc map for 
reclassification purposes, taking into account Table 4. 


The relative weights determinations of the route criterion, random index, Eigen values, and 


consistency of the three decision makers are shown in Table 5 — 7; while the summary of the weight 
selection by each decision maker is shown in Table 8. 
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Table 4: Reclassification of dataset 
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Table 5: Decision matrix weights and consistency of decision maker I 
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Table 7: Decision matrix weights and consistency of decision maker III 
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Table 8: Summary of weight selection by each decision maker 


Weight Selection (%) 





Table 9's Decision Maker I is a Principal Manager (Civil-Lines) who completed the decision matrix 
effectively, and when computed, the criteria slope of 25.49 percent is the highest relative weight and 
had the most influence on the HV-OHTL route selection. Decision Maker II, an environmental 
engineer, was tasked with identifying environmental impact assessment and recommending mitigation 
actions for the HV-OHTL project, and when the protected area of 26.57 percent was estimated, the 
decision matrix had a higher impact. The key worry of Decision Maker III, a maintenance engineer, 
was the availability of a good road network, which will provide access during the maintenance period 
with a 29.70 percent score. 


3.1. The cost distance and backlink raster layer 


A raster 1s initially created via a cost-distance analysis (Figure 9) that assesses the expense of 
travelling over a surface. The raster might utilize a single cost component, such as slope, or it could 
incorporate different sorts of charges. A cost direction raster, which advises a traveller which way 
incurs the least cost, and a cost accumulation raster, which keeps track of the overall effort expended 
as the traveller progresses, were developed using the cost raster. The least cost path between two 
points was calculated using these raster’s, showing the simplest route rather than the shortest. The cost 
distance and back link raster were used to do the least-cost path analysis. The source cells and the 
weighted sum are used to create the cost distance raster. The cost distance operation works in the 
background in a similar way to the Euclidean distance operation, except it calculates the shortest 
weighted distance between the cells. The back link raster specifies the pixel that should be used to 
return to the destination cell (ArcGIS Help 10.1). 
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Figure 9: The cost distance and backlink raster layer 
The route selected by Decision Maker I — III is represented in Figure 10, also the mean of the three 


decision makers was used to generate a fourth route. The total route lengths generated are shown in 
Table 8. 
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Figure 10: Route map showing preferred route of Decision Maker I-III and the mean route 
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Table 9: Length of routes based on Decision Makers 


Decision Makers Route Route Length Major Detours Max Average Slope 
(km) Slope (%) (%) 
Decision Maker I Route 1 174.79 30 11.1 0.9 
Decision Maker II Route 2 134.13 9 24.1 1.5 
Decision Maker II Route 3 174.36 32 10.3 1.1 
Mean of Decision Makers Route 4 138.05 20 10.0 1.1 


The number of angle points on each route, as well as the total distances, were taken into account while 
choosing the optimum route. The maximum and average slopes as they affect engineering design for 
each route were also taken into account. 


4.0. Conclusions 


Slope, Geology, Erosion, Soil, Land Use Land Cover, Road Network, and Protected Zone layers are 
part of the data set used in this investigation. To develop reclassified layers, these layers are put 
through a process. The applications for reclassification are carried out based on the viability of the 
layers' features. According to the HV-OHTL expert's understanding, the features are ranked by values 
ranging from | to 9. 


On HV-OHTL, there are three decision makers with various professional backgrounds. The AHP 
approach allows for the decision makers' preferences to be reflected. Because the route selection 
process for HV-OHTL is viewed from a distinct perspective, each decision maker prioritizes various 
layers. The cost surfaces are developed based on the preferences of each decision maker in order to 
describe the least expensive paths. The length, cost, impact aspects, and risk of the route possibilities 
are all considered. The cost of each route is calculated using a strong power line design program, 
power line system computer assisted and drafting, or PLS- CADD, that is generated from an actual 
HV-OHTL design. Material, installation, and shipping costs, as well as the cost of the cell where the 
power structure is located, are all included in the cost criterion. The number of intersections with 
linear structures is used to calculate the risk factor. The effect factor is calculated using the overall 
area (6707 km’) of the route right -of- way, which passes through settlements, agricultural areas, and 
forests. The decision makers’ consensus is used to calculate the weight of each criterion using the 
AHP pairwise comparison approach. The decision makers define the convenient preference functions 
and thresholds for each criterion in terms of its attributes. 


Recommendations: 


Prior to construction, the Power Transmission Company must be aware of all legal requirements for a 
new HV-OHTL. State laws and regulations must be observed and implemented to the letter. All 
stakeholders must be identified and their concerns addressed throughout the planning stage of a 
project. This will save time and money by avoiding unnecessary project delays. 


Because the resolution and precision of the data have a direct impact on the output information's 
reliability, the DEM layer should be generated using LiDAR data with a high resolution and 
operational flexibility. Small elevation changes will make the DEM more responsive. It 1s 
recommended that the numbers and background of decision makers be increased with a view to 
achieving unbiased and reliable solution for HV-OHTL routing decision problems. 


The model developed in this work should also be applicable to other voltage levels of power 
transmission lines route selection challenges, as well as other linear engineering structures. A national 
information centre or data bank should be formed with the goal of storing, updating, converting, and 
evaluating data that cabe used for research and other purposes. Civil-lines design engineer, 
environmental engineer, and maintenance engineer of HV-OHTLs projects are among the 
professionals chosen for consultation. They solely represent the Transmission Company of Nigeria's 
SMEs (TCN). Additional interviews with specialists from the Federal Ministry of Power, as well as 
officials from the Federal Ministry of Environment and landscape preservation organizations, would 
provide a more thorough assessment of the final route. 
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